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Abstract. Recent estimates of the properties of the Crab nebula are used to derive constraints on the moment of 
inertia, mass and radius of the pulsar. To this purpose, we employ an approximate formula combining these three 
parameters. Our "empirical formula" / ~ a{x)MR^, where x = [M / MQ){km/ R) , is based on numerical results 
obtained for thirty theoretical equations of state of dense matter. The functions a(x) for neutron stars and strange 
stars are qualitatively different. For neutron stars aj^s{x) = a;/(0.1 + 2x) for x < 0.1 (valid for M > 0.2 Mq) 
and a-Ms{x) = §(1 + 5a;) for x > 0.1. For strange stars ass (a;) = f (1 + a;) (not valid for strange stars with crust 
and M < 0.1 Mq). We obtain also an approximate expression for the maximum moment of inertia /max. 45 — 
(-0.37 + 7.12 • a:max)(Mmax/M0)(i?Af„„JlO km)^, where J45 = I/W'^ g • cm , valid for both neutron stars and 
strange stars. Applying our formulae to the evaluated values of /crab, we derive constraints on the mass and radius 
of the pulsar. A very conservative evaluation of the expanding nebula mass, M^eb = 2 Mq, yields Mcrab > 1.2 Mq 
and -Rcrab = 10 — 14 km. Setting the most recent evaluation ("central value") M^eb = 4.6 Mq rules out most of 
the existing equations of state, leaving only the stiffest ones: Mcrab > 1.9 Mq, Rcmb = 14 — 15 km. 
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1. Introduction 

The moment of inertia of neutron stars plays a crucial role 
in the models of radio pulsars. In the standard case of rigid 
rotation the total energy expenditure per unit time i^tot 
is related to the measured pulsar angular frequency Q and 
its time derivative by Etot — — /ri^i. However, even in 
the simplest case of slow rotation (il^ ^ [c/ Kf'GM / Rc^ , 
RVL <C c) characteristic of observed radio pulsars, the re- 
lation of / to the matter distribution within the star is 
complicated by the general relativistic effects such as the 
dragging of the local inertial frames. Among all global 
neutron star parameters the moment of inertia is the most 
sensitive to the dense matter equation of state (EOS). The 
"theoretical maximum mass" of neutron stars increases by 
a factor of two when going from the softest to the stiffest 
EOS. The "theoretical maximum moment of inertia" how- 
ever increases then by a factor of seven (Haensel 1990). In 
view of this particular sensitivity of / to the largely un- 
known equation of state (EOS) of dense matter at supra- 
nuclear density, it is of interest to look for observational 
evaluations /, which could then be used to constrain the- 
oretical models. 

In the present paper we use recent evaluations of the 
parameters of the expanding Crab nebula, in order to de- 



rive constraints on the Crab-pulsar neutron star. In Sect. 2 
we reconsider the energetics of the Crab nebula, and using 
expanding-shell model of the nebula we derive a constraint 
on the moment of inertia of the Crab pulsar. In Sect. 3 
we derive approximate formulae for I of neutron stars and 
strange stars, based on a statistical analysis of a numerical 
data sample obtained for thirty theoretical EOSs of dense 
matter. These formulae are then used in Sect. 4 to derive 
constraints in the M — R plane, implied by the value of 
-^Crab, deduced in Sect. 2. Appendix presents results of an 
analysis of the correlation between the maximum value 
of / and the mass and radius of the configuration with 
maximum allowable mass. 



2. Energetics of the Crab nebula and constraints 

on /crab 

The present values of and Cl for the Crab pulsar 
are 188.119 s^^ and -2.366 • 10"" s"^^ respectively 
(Lyne & Graham-Smith 1998). Therefore, the present rate 
of the rotational energy loss is 
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mVt = E,ot = 4.45 • 10^* /45 erg s 



(1) 
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It is commonly accepted that i^rot has to account for the 
luminosity of the pulsar and the nebula, as well as for the 
accelerated expansion of the nebula, so that 



ad 



E,- 



(2) 



Although the nebula has a complicated shape we will as- 
sume for simplicity that it is spherically symmetric with 
a radius i?ncb — 1-25 pc (see e.g. Douvion et al. 2001, 
Petersen 1998). The total power of radiation in the X, 
UV, and optical domains is given by Petersen (1998): 



E, 



rad 



1.5- lO-' 



2 kpc 



erg s 



(3) 



Nugent (1998) calculated the velocities of expansion of 
50 bright filaments by comparing their positions on four 
high-resolution photographs taken in 1939, 1960, 1976 and 
1992. After projecting the motion of filaments backward 
in time, he obtained convergence in AD 1130 ± 16 yr, 
in accordance with Trimble (1968). The calculated mo- 
ment of the supernova event is by A = 76 yr later than 
the date recorded by the Chinese astronomers, AD 1054. 
Therefore, the actual expansion had a non-zero accelera- 
tion. Assuming a constant value of the acceleration, v, and 
putting the present mean expansion velocity, determined 
from the measurements of the spectra, v 1.5-10^ cm s~^ 
(e.g. Davidson & Fesen 1985, Bietenholtz et al. 1991, 
SoUerman et al. 2000), we obtain 



2t;A 

rp2 



0.82 • 10~^ cm s"^. 



The Crab nebula expands in the interstellar, mostly 
hydrogenic, gas. Let the number density of the hydrogen 
atoms be nu- The expanding shell sweeps the interstellar 
medium, which becomes a part of the shell after being 
accelerated to the shell velocity. Consequently, the mass 
of the nebula increases at the rate 



(5) 



where mn is the mass of the hydrogen atom. The cor- 
responding rate of the kinetic energy injection needed 
to support the expansion is Mnebi'^/2. The total power 
needed to support the nebula expansion is therefore 



E,: 



(6) 



Under our assumptions, i?rot should exceed the sum 

-Brad 

the moment of inertia of the Crab pulsar 



i^exp- This condition results in a lower-bound on 



-?^Ciab > 



-^cxp ^" -^rad 



n\n\ 

which can be rewritten as 



-fCrab,45 > 0.34 



(7) 



2 kpc 



0.55 



neb 



0.15 
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neb 



Mq 
Tin 



1 pc y 0.2 cm~3 



(8) 



We will use the canonical distance to the Crab neb- 
ula given by Davidson & Fesen (1985), d = 1.83 kpc. 
The second term on the right-hand-side of Eq. (||) is the 
most important one - alas, it is also the most uncertain. 
Davidson & Fesen (1985) give the estimate 2 — 3 Mq, 
while MacAlpine & Uomoto (1991) evaluate this mass as 
1 — 2 Mq. The most recent estimate was obtained by Fesen 
et al. (1997), who get much higher value 4.6 ±1.8 Mq. 

In what follows, we will consider two estimates of Mncb- 
The lower one, 2 Mq, is a mean of the evaluations given 
by Davidson & Fesen (1985) and MacAlpine & Uomoto 
(1991). The higher one, 4.6 Mq, uses the "central value" 
of the mass calculated by Fesen et al. (1997). 

Most recent estimates of nn are based on the evalu- 
ation of the hydrogen column density in the direction of 
Crab. Such a method yields nn — 0.5 cm~^ (SoUerman et 
al. 2000, Willingale et al. 2001). However, this might be 
a severe overestimate of nn, because the Crab radiation 
goes through the Galaxy disk on his way to the detector. 
In view of this, we will prefer to use a conservative value 
TiH = 0.2 cm~3 quoted by Manchester & Taylor (1977). 

Using the above mentioned values of the quantities en- 
tering the right-hand-side of Eq. (||), we get two lower 
bounds for the moment of inertia /crab, corresponding 
to two different estimations of the Crab nebula mass: 

-^Crab,45 > 1-61 for Mncb = 2 Mq, and /crab,45 > 3.04 



for Mn 



4.6 M, 



0- 



(4) 3. The I{M, R) correlation 



In their seminal paper, Ravenhall & Pethick (1994) stud- 
ied the properties of the ratio I/MR^^, where i?oo stands 
for an apparent (radiation) stellar radius, determined in 
the analyzes of thermal radiation from the close by neu- 
tron stars. 



Rno = 



R 



^\ - 2GMlRc^ 



(9) 



where R is the coordinate radius of the space-time metric. 
In particular, they found that for an intermediate range 
of (M/MQ)(km/i?) - 0.05 - 0.15 the ratio //Mi?^ stays 
nearly constant, and depends rather weakly on the dense 
matter EOS. This finding resulted in their approximate 
formula for /. 

In our quest for a universal, i.e., independent of the 
specific EOS approximate relation between /, M and 
i?, we relied on a large set of more than thirty EOSs 
of dense matter. This set includes various models of 
dense matter, which were obtained under different hy- 
potheses on the composition of the matter at p > 2pQ, 
where po = 2.7 • lO^'' g cm^^ is the normal nuclear den- 
sity. Considered compositions of the super-dense mat- 
ter in neutron stars (NS) include: nucleons and leptons 
(EOSs of Baldo et al. 1997, Bombaci 1995, Balberg & 
Gal 1997; Balberg et al. 1999; Bethe & Johnson 1974, 
Pandharipande 1971, Pandharipande & Ravenhall 1989, 
Douchin & Haensel 2001, Wiringa et al. 1988, Walecka 
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1974, Haensel & Proszynski 1980, Akmal et al. 1998); nu- 
cleons, hyperons, and leptons (EOSs of Glendenning 1985, 
1997, Balberg & Gal 1997, Weber et al. 1991, Balberg et 
al. 1999); deconfined quark matter mixed with baryonic 
matter (EOSs of Glendenning 1997); pion and kaon con- 
densates coexisting with baryonic matter, in a pure or 
a mixed-phase configuration (EOSs of Glendenning 1997, 
Kubis 2001). We considered also three models of self- 
bound absolutely stable strange quark matter, forming hy- 
pothetical strange stars (SS; EOS of Haensel et al. 1986). 
After studying various types of correlations between /, M 
and R, we found, in accordance with Lattimer & Prakash 
(2001), that it is the ratio I/MR^ which has the most uni- 
versal dependence on the compactness M/R. Our result 
is seen in Fig. ^ 




M/R [Mg/km] 



Fig. 1. The ratio I/MR^ versus compactness parameter 
M/R for selected set of EOSs of dense matter. Branches 
NS and SS correspond to neutron stars and strange stars. 
Thick lines were obtained using the analytical fitting for- 
mulae (solid - NS, dashed - SS), described in the text. 
Both thick lines are terminated at the M/R upper-bound 
for the EOS respecting condition t^sound < c, which is equal 
0.24. 

In what follows, we will use, following Ravcnhall & 
Pethick (1994), the dimensionless compactness parameter 
X = M/R-km/MQ. The dependence oil/MR^ on x is well 
approximated, up to M^^^, by a function a{x), different 
for neutron stars (NS) and strange stars (SS), 



aNs(a;) = 
ass (a;) = 




The difference in the x-dependence of I /MR? between 
neutron stars and strange stars stems from the qualita- 
tive difference in their density profiles. The density pro- 
file within a moderate mass strange star is very flat. On 
the contrary, density falls by many orders in the outer 
layer of a neutron star. Consequently, at a same x < 0.2, 
(//Mi?2)ss > iI/MR^)m (the case of a; > 0.2 will be 
discussed later in this section). As we see in Fig. ^ the 
precision of the formulae given in Eqs. (^ is typically bet- 
ter than 10%. The formula for neutron stars breaks down 
only at very low M, below 0.2 Mq. Such low masses do not 
seem relevant for neutron stars observed as pulsars. |^ The 
formula for ass looses its validity only for very low-mass 
{M < 0.1 Mq) strange stars with crust. On the contrary, 
the formula for / becomes exact for low-mass bare strange 
stars, because they are Newtonian and their density is to 
a very good approximation constant, which corresponds 
to ass - 2/5. 

The overall upper-bound on x is reached at Af,„ax for 
a pure causality-limit EOS P = c^{p — ps), and does not 
depend on p^: a^^ax = 0-240 (see Lindblom 1984, Haensel 
et al. 1999). Within General Relativity, static compact 
star models with EOS satisfying iisound < c have x < aj^^x- 
We obtain aNs(2:Sax) = 0-489 and ass(a;max) = 0.496, so 
that in the x — x^ax limit ass — aNS within about 1%. 

4. Estimate of /crab aid constraints in the M — R 
plane 

In Sect. 2 we obtained two lower bounds for the Crab pul- 
sar moment of inertia /crab , which depend on the different 
evaluations of the expanding shell mass: /crab,45 > 1-61 for 
Mnob = 2 Mq (conservative estimate), and /crab,45 > 3.04 
for Mnob = 4.6 Mq (newest estimate). 

These bounds on /crabj combined with our fitting for- 
mulae, derived in Sect. 3, imply constraints on the mass 
and radius of the Crab pulsar (for details, sec Fig. ||). 
Neutron star configurations below the thick solid lines are 
ruled out. The weaker, conservative constraint eliminates 
a number of soft and moderately stiff EOS. 

For those EOS which are acceptable, this constraint 
gives the range of possible mass and radius of the Crab 
pulsar, Mcrab > 1-2 Mq and i?crab = 10-14 km. 
The second (Mneb = 4.6 Mq) constraint is much more 
dramatic. Only a very stiff EOSs are accepted, because 

Mcrab = 1.9 Mq and i?Crab = 14-15 km. 

For the strange-star model we obtain strong con- 
straints (thick dashed line in Fig. 2) even for the conser- 
vative estimate Mncb = 2 Mq: m'^^^X > 1-5 ^'^0)-Rcrab = 
10—11 km. Putting A/nob — 4.6 Mq excludes definitely the 

^ Approximate formula for very low values of M/R corre- 
sponding to the vicinity of the Mmin — 0.1 Mq is given in 
the footnote on p. 850 of Ravenhall & Pethick (1994). (We be- 
lieve, that we have found a couple of misprints in this foot- 
note: the equation for very small x — M/R should be rewrit- 
ten as 7 ~ 7.5 • M^R, and for greater x, I/MR^K[R) ~ 
0.23 - 0.07 • (1 - lOxf ■ V^) 
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R [km] 



Fig. 2. The mass-radius curves for a large sample of EOS 
(thin dotted curves - neutron stars, thin dashed curves 
- strange stars) and constraints resulting from the esti- 
mates of / of the Crab pulsar (thick solid line - neutron 
stars, thick dashed line - strange stars), obtained in the 
text. Hatched area is excluded by the General Relativity 
combined with Usound < c. 



strange-star model. Note, however, that the strange-star 
constraint is likely to be an academic one anyway because 
the glitches in the timing of the Crab pulsar are widely 
believed to rule-out the strange-star model (Alpar 1987, 
nevertheless see Glendenning & Weber 1992). 

The formula expressing / in terms of M and R may 
be useful for putting constraints, resulting for pulsar ob- 
servations, on the neutron-star models. Consider a pulsar- 
powered supernova remnant of the Crab type. If the mea- 
sured total power needed to sustain the radiation of the 
system and the acceleration of the nebula expansion is 
i^totj this implies the lower bound on the pulsar moment 
of inertia, / > i?tot/(f^|f^|), where 17 and are measured 
pulsar angular frequency and its time derivative. Reliable 
estimates of -Etot, combined with "empirical formula" for 
/, would then imply constraints on the neutron-star mod- 
els in the M — R plane. Bounds for the mass and radius 
of the Crab pulsar obtained for M^ch — 4.6 Mq are very 
strong and should therefore be considered with caution. 
More advanced models of the structure and evolution of 
the Crab nebula deserve to be developed. This problem 
will be addressed in a subsequent paper. 
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Appendix: Correlation between /max and M^aax, 

The fact that cns and ass converge at the highest values 
of compactness parameter suggests that there might ex- 
ist a universal correlation between the maximum value 
of /, /max, and Minax and i?M„ax> valid for both neu- 
tron stars and strange stars. Such a correlation has been 
pointed out by Haensel (1990) for a set of nineteen 
EOS developed in 1970s and 1980s: /max,45 = 0.98 • 
(Mmax/M0)(i?M„ax/lO km)2. For the EOS set consid- 
ered by Haensel (1990) this "empirical formula" repro- 
duced values of /max with an overall accuracy = 
0.03. Within our sample of numerical results, obtained 
for selected thirty EOSs, we found correlation /max — 
C(Mmax/M0)(i?M„,a^/lO km)^, the same as in (Haensel 
1990), but with a slightly lower numerical factor C = 
0.97; the overall fit is characterized by = 0.05. The 
value of can be significantly reduced if one takes 
into account an additional correlation between the ratio 
/max/(A/max/?M^^^) and the maximum compactness pa- 
rameter ajmax reached at Mmax- The improved empirical 
formula yields: 

/max.. ^ (-0.368+7.122..max) (^) {^)^ .(H) 

It is valid both for neutron stars and strange stars, with 
the overall quality of the fit x^ ~ 0.02. The coefficient 
C = 0.97 in the simplest empirical formula corresponds to 
an "average value of Xmax", aimax = 0.19. Some general 
statements concerning different classes of the EOS can be 
made. The "empirical formula" , Eq. (|ll|) , under-estimates 
/max for the superluminal EOSs. On the contrary, it over- 
estimates /max for the EOSs with high-density softening 
due to the appearance of hyperons. Finally, for the "mini- 
mal composition EOSs" (nucleons and leptons) as well as 
for the strange matter EOSs, our formula reproduces /max 
within better than 5%. 

A different, very precise formula for /max of strange 
stars (bare or with crust) results from "scaling proper- 
ties" (Haensel et al. 1986) of the configurations with Mmax 
and /max- For the simplest model of strange quark mat- 
ter composed of massless, non-interacting quarks, scaling 
with respect to the value of the bag constant implies ex- 
act relation /max,45 = 0.94 (Mmax/M©) (i?Af„,J10 km)^ 
After inclusion of the QCD interaction and finite strange- 
quark mass the above equation is no longer exact. For 
rrisC^ = (150 — 200) MeV, the numerical prefactor to be 
put in the expression for /max, 45 is 0.97, and the formula 
is precise within a few percent. 

One might argue that the validity of a simplest "em- 
pirical formulae" for /max is a consequence of the fact 
that realistic EOSs are well approximated by the poly- 
tropes. For a polytropic EOS P = Kn\^ (where rib is 

^ In the present paper the value of is actually scaled by 
the variance, i.e., it is equal to the standard value of divided 
by . Such a treatment is convenient in our case, where we 
have to compare the qualities of different fits. 
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baryon density of matter) with a fixed F, the relation 
-fmax,45 = C ■ {M/Mq){Rm^,^/W km)^ is exact, with C 
independent of K. We get Cr=2 = 0.83 and Cr=a = 1.04, 
so that the "best fit" value for realistic EOSs lies roughly 
in the middle of the F = 3 and F = 2 values. However, 
the realistic EOSs at supcr-miclcar densities are not poly- 
tropes. One can speak only about a local, density depen- 
dent F(nb). For example, in the case of dense matter with 
hypcrons the EOS softens considerably (and in a discon- 
tinuous way) at the thresholds for the appearance of a 
new hyperon, with F dropping below one (Balberg & Gal 
1997). Clearly, the validity of the empirical formula for 
-''max is not due to the polytropic character of the EOS of 
dense matter. 
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